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Abstract 

It is known that any minimal model Mp receives along its (psi irrelevant di- 
rection two massless integrable flows: one from Mp+i perturbed by 4>i3, the other 
from Zp-i parafermionic model perturbed by its generating parafermion field. 
By comparing Thermodynamic Bethe Ansatz data and "predictions" of infrared 
Conformal Perturbation Theory we show that these two flows are received by 
Mp with opposite coupling constants of the (p^i irrelevant perturbation. Some 
comments on the massless S-matrices of these two flows are added. 
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1 Introduction 



Many two-dimensional massless Quantum Field Theories can be seen as Renormaliza- 
tion Group (RG) flows connecting an ultraviolet (UV) Conformal Field Theory (CFT) 
to an infrared (IR) one When an infinite number of charges is conserved the theory 
is integrable. In this case a powerful method of investigation consists in the Thermo- 
dynamic Bethe Ansatz (TBA) 0, that can be thought as a set of non-linear coupled 
integrable equations driving the RG evolution of various quantities along the flow. In 
principle the TBA equations could be deduced from the (massive or massless) S-matrix 
of the theory, but the derivation is technically easy only in the case of diagonal S- 
matrices. If a kink structure connecting different colored vacua appears, one must 
resort to very complicated Bethe Ansatz techniques to diagonalize the color transfer 
matrix and deduce the TBA system. In spite of this, a lot of TBA systems attached to 
many integrable theories have been conjectured, and extensively checked against UV 
perturbative results. In the case of massless theories, however, also the IR limit is given 
by a non-trivial CFT to which the RG flows is attracted by an irrelevant integrable 
operator $(a;). This simplified sentence means that one can write an IR-effective action 
approximating the behaviour of the theory at large scales 



where Air (formally) represents the action of the IR CFT. The "h.o.t." means higher 
order terms, i.e. an infinite number of higher dimension operators that in principle can 
contribute as counterterms each one with its independent coupling constant. In other 
words, this effective theory is not renormalizable. This makes the IR perturbation the- 
ory much more difficult than the UV one. In these two-dimensional theories, experience 
tells us that the UV perturbative series usually has a finite radius of convergence. In- 
stead, the IR one is at best an asymptotic series. In spite of these difficulties we shall 
see in this letter that the first few orders of the IR series can be reasonably controlled 
and turn out to be of interest when compared with the exact results coming from the 
integration of the TBA equations. This allows further checks of the validity of TBA as 
well as calculation of quantities that can shed new insight on the structure of the space 
of RG fiows, as we shall see below. 

The idea of comparing TBA results with IR perturbation theory traces back to 
Al.Zamolodchikov |^, but only the case of TT perturbations (T being the stress- 
energy tensor) were developed enough. Klassen and Melzer went a bit further, 
by exploring the simplest cases of 03i IR perturbations of minimal models. Much of 
the inspiration of the present work comes from their approach. IR perturbations have 
been dealt with in a more systematic way, although in a rather different problem, by 
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Berkovich Also this paper contains elements that have been illuminating for our 
analysis. 

In this letter we make use of IR perturbation theory to compare the TBA results of 
two celebrated examples of integrable massless field theories, maybe the most studied 
ones: 

1. the minimal models Mp+i perturbed by 0i3 that notoriously flow to Mp. Follow- 
ing 1^ we denote these massless theories by A^Ap^\. 

2. the Zp_i-parafermionic theories perturbed by their = ipiipi ipl^pl operator {ipi 
being the generating parafermion) , which also are known to flow down to Mp. 
Following they will be denoted in this paper by 

That the minimal models and the parafermion theories have very strict relation in gen- 
eral, even at the pure conformal level, is not a new surprise and can be traced back to 
the fact that both heavily involve the Ai'^ affine algebra in the deep of their construc- 
tions. The interesting fact here is that both theories flow down to Mp "attracted" by 
the same operator 03i. One can wonder which is the feature of the IR 03i perturbation 
theory that distinguishes the two flows. Our investigation has the aim of clarifying this 
issue and add a new piece of information to the beautiful puzzle of understanding the 
map of 2 dimensional integrable flows. 



2 TBA equations for MA^^^^ and i/^^) models 

A set of TBA equations for AiA^^ has been proposed by Al.Zamolodchikov Q], and 
checked against various tests that give more than reasonable confidence to its cor- 
rectness. Basically the prediction for the scaling central charge c(r) (r = MR is an 
adimensional quantity parametrizing the RG flow in terms of the mass scale M and 
inverse "temperature" R) is 

Q n+OO 

c{r) = - Y^dOv^Ue) (2) 

71 J — oo (J 

where Uai^) = ^^^^a,i + ^e~^6a,p-2, the index a running from 1 to p — 2 along a Ap_2 
Dynkin diagram. The La{9) = log(l + e"^"*^^-*) are evaluated as solutions of the set of 
integral equations 

1 f + OO 

ya{e) = ea{9) + — / d9'v{e - 9') Y LbU{9') (3) 
Zn J-oo 

with kernel if{9) = 1/ cosh(^^) and lab incidence matrix of the Ap^2 diagram. With 
this TBA one can test even for low p the validity of the results on the AiA^^ models 
established perturbatively by A.Zamolodchikov for high p. 
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A similar TBA set has been proposed by Fateev and Al. Zamolodchikov 0] to 
describe, or even better, to give the first evidence, of highly nonperturbative flows 
between parafermions and minimal models, namely the Hj^^^ theories briefly described 
above. Again, the TBA equations and the formula for the scaling central charge are 
given by eqs.(|D and @, with the same kernel, the only difference now consisting in a 
different lab incidence matrix, namely that of the Dp Dynkin diagram. 

3 IR perturbation theory 

The discussion of IR perturbative series is basically plagued by the existence of infinitely 
many irrelevant operators Oj{x) that can contribute to it each one with its independent 
coupling Qj, thus making the theory not renormalizable. However, conformal perturba- 
tion theory has the advantage to be strongly constrained by the fact that the operators 
appearing in the correlators in the perturbative series must be found among the oper- 
ator content of the unperturbed CFT itself. The effective perturbed theory is defined 
by the action 




We use the following requirements to restrict the space of IR perturbing operators V 

• all perturbing operators must be scalars, i.e. Aj = Aj in order to keep 1+1 
Lorentz invariance 

• all perturbing operators must be irrelevant, Aj > 1 

• total derivative operators do not contribute to the action 

Even after these restrictions, the space V is still extremely large and difficult to analyze. 
But there is another issue which is the main characteristic of these perturbations: their 
integrability. In other words, we know that these perturbations describe, near IR, flows 
having an infinite number of conserved currents. In particular it is known that both 
series of models AiA^^ and H^^'^ have for all p, an infinite set of local integrals of 
motions (IM) of spin 1 mod 2. This does not exclude for some specific models in these 
series existence of other IM's, but at least this series of IM's must be kept conserved 
all along the flow. All operators destroying even only one of these currents must be 
discarded from V. This is the really powerful constraint we can put on our series. It is 
indeed a fact stressed by many authors that only operators from the families of [1] , [(pis] 
and [4>3i] have the nice property to preserve all the IM's of spins 1 mod 2. Therefore 
our space V can be restricted to these 3 families in the following analysis. 
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If R denotes the radius of the cyhnder where the theory is put for TBA analysis, 
the scahng function (normahzed as to give the UV and IR central charges in the two 
limits i? — i> and R oo respectively) is given near IR by a series expansion in the 
correlation functions of the Oj{x) operators on the cylinder 

oo 

n=l ji,...,j„6V 

where 

Pn,-,. = l2^—^R'-^-^y^^ / l[d'^^^(0,MI[0,A^^^))cyl (6) 

i=2 i=2 

The coupling constants Qi carry a scale dimension ?/i = 2(1 — Aj) each, and as in the 
usual UV case, there will be a relation with the mass scale M of the form 

= n,My^ (7) 

where Ki are numerical coefficients. Comparison with TBA data can be done by first or- 
dering the Z^iLi Ujt — Ua in decreasing order (all yi are negative for irrelevant operators). 
Define 

— E -^iivJn'^ji---'^jn (8) 

Ja 

where the set J a is formed by all the possible choices of ji, jn such that Z^ILi Vh ~ Va.- 
then the scaling function c(r) has the expansion 

c(r) = ciR + Yl Car^" (9) 

a 

that can be directly compared with the numerical integration of the TBA equations. 
What in practice one has to do is to choose a truncation of this series and try to fit it 
against the TBA numerical data. 

In order to explicitly compute the Pjj^,...,j„ one has to transform them from the 
cylinder to the sphere. Note that in the case of secondary operators the transformation 
can be more complicate than the usual one for primary operators considered in [§] to 
give the "sphere" formulae for P's. Luckily, as we shall see, for the purposes of this 
paper, the only secondary operator we have to consider is the TT secondary of the 
identity. The well-known transformation of the stress-energy tensor creates a term ^ 
that allows a contribution from the 1-point function of this operator on the cylinder. 
One must keep in mind this exception in the following analysis. No secondary operator 
of families different from [1] can have a non-zero cylinder 1-point function. 

In writing down the asymptotic series (y) one has to analyze which correlation 
functions really give a contribution and to which order. Here we give a brief summary 
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of this analysis in the case of the three famihes of integrabihty preserving operators 
mentioned above. 

The [1] family contains as first (scalar) secondary the determinant of the stress- 
energy tensor TT. Its renormahzation group eigenvalue y = 2(1 — A) is —2, therefore 
we expect contribution from its 1-point function at in the IR series. Its 2-point 
function will contribute at r~^, etc... Next operator in the identity family is : T^T^ . 
whose first contribution appears at order r~^. 

The operator is the dominating one in the IR series. The first contribution 
comes from its 2-point function (03i03i) and, as A31 = 1 + ^, which implies y = —4:/p, 
it is at order r~^^^. Next contribution is given by the 3-point function at r~^^/^. This 
will be the crucial point of our analysis, as we shall see. Further n-point functions of 
031 appear at orders r^'^^lv. Secondaries of 03i can give contributions at orders r~^~^/^ 
or higher, as it appears from analysis of the "lowest" case, namely (03i03i )• (Here 03i 
denotes a generic secondary of 03i at level 2). The 03i expansion is unaffected by its 
secondaries if we do not consider terms 0(r^^). Another interference that can appear 
is between 03i and [1] family. The first term comes from (03i03iTT), this is at order 
and if we consider only terms of order less that this we are also free of this 

problem. 

Let us now turn to the contributions from the pure [013] family. As ^13 is relevant, 
the dominant field of this family allowed to be in V is L_2Z_20i3. Its 2-point function 
contributes at order 7^-6+8/ jj^ g^^y ^^^q higher than the limit we already imposed to 
ehminate interference between ^31 and TT. All other correlators of secondaries of ^13 
contribute at order higher than this. Insertions of secondaries of the identity into such 
correlators raise the order even more. Finally the interferences between secondaries of 
013 and 031 also happen at higher orders, because the 2 and 3-point functions involving 
both 013 and 03i are all zero. 

We conclude that, if we are content to consider contributions up to r~^~^/^ excluded, 
we can only consider the two operators 03i and TT in the perturbation expansion. 
Higher orders not only become more complicate from the "theoretical" point of view, 
but a careful numerical analysis shows that their contributions begin, for a typical 
numerical sample of TEA data in the range 200 < r < 2000, to be in the region where 
the dehcate equihbrium of asymptotical convergence starts to break. Therefore we 
consider in the following the effective action 

A^Am + j (fx{ghi + tTT) (10) 

To finish this section, we have to mention another delicate issue that appears in the 
perturbative series. When some operator is such that for some integer it happens 
that Ny = —2, the corresponding A^-point function is divergent, even after analytic 
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0.04730(6)r-9 + 0.194(3)r-3 + 0.98(6)r-^ - 0.0(l)r-2 - 1.2(5)r-^ 



Table 1: The numerical fits of TBA data against IR perturbative series 



continuation in y. To avoid this singularity, one should introduce a counterterm that 
results in the appearance of a logarithmic term of the form r^^ logr that must be taken 
into account when fitting with numerical data. This happens in the model Mp with p 
even for the |-point function of 03i. Therefore, in the M4 model this problem arises 
already at the first term, i.e. the 2-point function of ^31. Similarly, in the Mq case the 
3-point function (which is the crucial one in our analysis) is plagued by this log-term. 
However, for higher Mp with even p, this interference moves to higher and higher ^31 
correlators, and leaves the 2 and 3 point functions free of this problem, allowing clean 
application of our argument. 

4 Numerical analysis 

We have performed a high precision (14 significant digits) numerical integration of the 
TBA equations in the range 200 < r < 2000 step 10. Stability of the results has been 
checked against increasing the number of iterations and decreasing the integration step 
in 9. The data in this range have been fitted with a function of at the powers 
predicted by the IR perturbation analysis done in the previous section, truncated at 
around r~^-^. The reason for this choice is explained a httle below. The coefficient of 
the fitted expansion thus found are reported in Tab.l. 

We started our analysis with the M5 IR theory. The M4 case presents the mentioned 
anomaly to have the first coefficient divergent, which results in a logarithmic contribu- 
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0.037121126 


0.194123 



Table 2: The theoretical values of the coefficient determined by the 3-point function of 
</)3i. The values of k needed for this calculation are deduced from the coefficient of the 
2-point function term. 

tion. On the other hand, it has been extensively discussed by Klassen and Melzer in . 
They give arguments to support the hypothesis that the two flows from M5 and Z3 come 
into M4 from opposite values of the 03i coupling g. Log-terms are totally absent in the 
next model M5. The dominating coefficient comes from the 2-point function of 031. It 
can be evaluated theoretically by resorting, e.g. to formula (2.28) of 0, which is valid 
also for the irrelevant 03i because it is a primary field. Out of that formula we extract 
the coefficient n relating the 0i3 coupling g to the mass scale, \g\ = kM~^^'p. Then use 
of this K and of the structure constant C*(|i j (31) that one can get from [T^ allows to 
compute the "theoretical" value of the coefficient of the 3-point function. The only case 
where this does not work is Mq where the 3-point function of 03i diverges leading to 
appearance of a log-term and also interferes with TT. In all other cases we are able to 
report the theoretical coefficient that we use as a check of the numerical correctness of 
our data. The most important thing that one immediately sees from Tab.l is that the 
coefficient of the 3-point function of 03i in the AiA^i model is always opposite to that 
of theif^' l\ one, a clear sign that the two flows enter the IR model from opposite 
directions! This definitely settles the problem of the two 03i attracted flows that enter 
each minimal model. But a more careful analysis of the numerical data of Tab.l show 
other interesting surprises. For example, within numerical error it seems plausible to 
take, at least in all models with odd Mp as IR, the coupling with TT equal to zero for 
the parafermionic flows, while it is different from zero for the minimal models ones. If 
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this fact has a relation with the different forms (with of without 9 shift - see below) of 
the two massless S-matrices of the models is an open and interesting issue. For even 
p instead, it happens that the TT coupling always interferes with some 'H.-point 03 1 
function. It is curious to note that the renormalized couplings still show the desired 
behaviour (equal absolute value and opposite/equal signs for odd/even n). 

Our numerical data seem to be reliable up to some r~^'^. Introduction of terms 
higher than this do not seem to improve the fit, in some cases even they make conver- 
gence worse. This is a typical behaviour of asymptotic series, adding new terms can 
give worse approximations of the function after a certain point. It could be possible 
that this r~^'^ plays a role of a limit after which asymptotic convergence becomes more 
problematic. Luckily, this also is the limit above which secondaries, interferences be- 
tween different families and other such phenomena are confined (at least for p < 16), 
a lucky situation indeed allowing us to observe the sign change in g so important from 
the physical point of view, before it becomes hidden by all these more or less uncon- 
trollable contributions. Note that number of 03i correlators which fall before this limit 
increases with p. Curiously, higher minimal models allow a better IR analysis than the 
lower ones. Indeed starting from My we are able to see also the 4 and 5-point functions. 
The 4-point function, as expected, gives the same coefficient in the two flows, while the 
5-point again presents the change of sign, thus confirming once more that the two flows 
come to Mp from opposite directions. 



5 Some comments on the S-matrices 

For the two classes of flows under considerations, massless S-matrices can be written. 



The minimal models case has been analyzed in |ri|, in relation with the so called 
"imaginary Sine-Gordon" model of relevance for polymer physics. The Sll = Srr 
matrix of Aip is formally equal to the S-matrix of the Mp_i model perturbed by 0i3 in 



the massive direction, and therefore it is given by the one described in |T^. The Sjil 



and Sir blocks, instead, turn out to be proportional to the Sll block, calculated for a 
value of 9 shifted by an imaginary quantity a. 

SLR{9)(xSLL{9 + ta) , SnL{9)(xSLL{9-ta) (11) 

The shift can be fixed by requiring that the two "mixed" components are equal, which 
should correspond to requiring time-reversal invariance. This imposes « = (a* = 
^^), n = 0, 1, .... As shown in ref. |11] the correct choice for this case is given by a = 



The proportionality factor is fixed by the analytical crossing- unitarity requirement. 

For what concerns the H^'^^ flow, the massless S matrix can be obtained by the 
following argument. Fateev wrote the massive S matrix of the H^^^ theories, i.e. the 
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Zp-parafermionic models perturbed by their fundamental parafermion operator in the 
massive regime. He got the S matrix as quantum group reduction of the S matrix of the 
Sausage model (see ref. FOZ). This latter is basically the spectral parameter dependent 
R-matrix of the sl{2)q quantum group in the spin 1 representation (in contrast to the 
Sine-Gordon S matrix, which is in the spin ^ representation). At q root of 1, the 
quantum group reduction gives the S matrix of the H^^^ models, as one can easily check 
against the known cases p = 2,3, 4. If this procedure works fine in the massive sector, 
we do not see any reason preventing its extension to the massless regime. The S-matrix 
for the massless Sausage model is written in [Q, and consists of blocks all equal to the 
Sine-Gordon S-matrix, i.e. Sll = Srr = Srl = Slr = Ssg- This is not unexpected, 
if we remember that the IR limit of the massless sausage is the c = 1 CFT, that is 
also the UV limit of Sine-Gordon theory. The g-reductions of this S-matrix describe 
the parafermionic flows H^^\ This also implies that the Sll and Srr components are 
equal for the two flows AiA^i and -ffp!^\. Again, this is not surprising, as the two 
flows share the same IR limit Mp. The distinction between the two cases is given by 
the Slr = Srl blocks that in the A4A^p^i case have the shift in 6, eq.(|ll|) while in the 
case SLRie) = SLL{e). 
This observation would exhaust the problem of writing the massless S-matrices for 
the two series of flows examined in this paper. However, one has to check the coinci- 
dences H^''^=MAi^^ (Ising model) and i7^^^=A^A^+^ (3 state Potts Model). While the 
first one does not give any problem, the Potts case is rather intriguing. The Sll and 
Srr blocks coincide, as expected, but the Slr = Srl do not. On the other hand, a TEA 
analysis similar to the one done in [|ll| (and generalizing the massive case of for the 
scaling tricritical Ising model) shows that both S matrices give the correct cuv = I and 



CiR = ^ (note that that the S matrix with unshifted Slr giving cuv = I in [|TT] differs 



from ours by a non-trivial Z-factor). It is evident that there must be a transformation 
of the asymptotic particle basis relating one S matrix to the other. The most appealing 
feature in this respect is the fact that the two theories AiA^^^ and -^3'^'' are in a sense re- 
lated one another by an orbifold procedure. Consider their UV limits: the AiA^^^ is the 
diagonal minimal model M5, the -^3'^'* one is the Zs-parafermion, which coincides with 
the complementary modular invariant {{A4,D4) in the Cappelli, Itzykson and Zuber 
classification ||l5l). One passes from one to the other modular invariant by an orbifold 
procedure. If we make the hypothesis that also off-criticality the two models are related 

(it) 

one another by such an orbifold, the S matrix of the model will be obtained from 
the AiAl^^ one by the orbifold recipe explained in This indeed works fine in the 
massive analogs of these two models, as well as for the Sll = Srr blocks here, which 
turn out to be equal in the two cases, thanks to the automorphism A3 = of their 
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adjacency diagrams. However, straightforward application of this folding procedure to 
the Slr = SfiL blocks would predict also for them equality in the two models, while, as 
seen, they are definitely not equal in this case. Our opinion is that the folding procedure 
a la Fendley Ginsparg [|1^ has to be modified slightly in the massless case, in the sense 
that the role of left and right movers are inverted: if we denote by — , 0, + the 3 nodes of 
the i?+ right mover behaves exactly as L_ left mover in the folding procedure, and 
so on. The — ^ + symmetry of ^3 guarantees that if only right or only left objects 
are concerned {Sll and Srr) this inversion is uninfiuent. If instead we are considering 
the scattering of a left and a right movers this modifies the folding procedure in such 
a way to give the shifted Slji from the unshifted one (they are related one another by 
sending sinh i cosh in this case. We intend to return to this curiosity and to the 
general problem of writing S matrices for non-diagonal 0i3 perturbed minimal models 
in massive and massless regimes in future. 



6 Conclusions 

In this letter we have given evidence that the two flows entering any minimal model 
from 031, namely the integrable theories AiA^^ and H^\, come into their IR limit 
from opposite directions. This result can be interesting in itself, as it sheds further 
light on the map of the two-dimensional RG space of integrable flows and of minimal 
models in particular. Further things should be investigated to better understand this 
issue. We have seen in the previous section that the "fine structure" distinction between 
different modular invariants may play a role in solving apparent contradictions between 
S matrices that appear to be different and are expected to be equal, as in the commented 
h'^'' versus M.A^^^ case. Also, the flows examined in this paper are not unique: the 
parafermions themselves have in general more than one modular invariant. It would be 
interesting to solve the whole puzzle of the flows between all parafermions and minimal 



models ones. This could also point out, by "undoing the truncation" procedures |TT|, to 
some orbifold version of the imaginary Sine-Gordon and of the Sausage models. Work 
in this direction is in progress. 

We think that the most appealing fact of the result of the present paper is the way it 
has been obtained, i.e. by use of IR perturbation, an instrument used only marginally 
so far in TEA analysis. It is possible that a better development of this IR techniques, 
maybe by comparison with the similar ones developed in P|, could provide new tools 
of investigation of the massless integrable theories. 
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